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O Abstract 



We report the results of an experiment to determine whether the half-life of '^'^Au depends on the shape of the source. This study 
was motivated by recent suggestions that nuclear decay rates may be affected by solar activity, perhaps arising from solar neutrinos. 
If this were the case then the y6-decay rates, or half-lives, of a thin foil sample and a spherical sample of gold of the same mass 
^ and activity could be different. We find for '^**Au, (ri/2)foii/(ri/2)sphere - 0.999 + 0.002, where T1/2 is the mean half-life. The 
maximum neutrino flux at the sample in our experiments was several times greater than the flux of solar neutrinos at the surface of 
the Earth. We show that this increase in flux leads to a significant improvement in the limits that can be inferred on a possible solar 
contribution to nuclear decays. 
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I Extensive research has shown that the decay rate of a radioac- 
tive nucleus is independent of its environment, except in those 
instances involving electron capture, internal conversion, or 
high external magnetic fieldsjilll. Therefore, it is surprising that 
a recent series of papers have raised the question 

of whether nuclear decay rates as observed in terrestrial labo- 
ratories are being influenced by solar activity. These authors 
point to correlated seasonal variations of the decay rates of ^^Si 
at Brookhaven National Laboratory (BNL) [7], and of '''^Eu, 
'^"^Eu, and ^^^Ra at the Physikalisch-Technische-Bundesanstalt 
(PTB) in Germany |8, 9]. In addition, they found a statistically 
significant drop in the counting rate of ^"'Mn coincident in time 
with the solar flare of 2006 December 13 HHIl]. 

Although it may be the case that the observed fluctuations in 
the BNL/PTB data reflect some poorly understood instrumen- 
tal effect, such as a seasonal dependence of detector efficiencies 
on external temperature, in the absence of additional evidence 
it remains an experimental question as to whether or not solar 
activity is in fact affecting the intrinsic decay rates. If this were 
the case then we might expect a range of responses from differ- 
ent nuclides arising from the same details of nuclear structure 
which account for the wide range of observed half-life {T1/2) 
values. Given the unknown origin of the effect, it is thuspossi- 



ble that analyses of other data motivated by Refs. Hll, 311,11] 
have either not seen 

lEintl orobs 
cay rates consistent with this work. 



by Ke 

have either not seen iflolfTlll or observed 1I12I1 time- varying de- 
To account for the data in Refs. ||2> B B IH Hi, it was pro- 
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posed that unstable terrestrial nuclei are being influenced by an 
(as yet unknown) interaction with solar neutrinos Ja]. A di- 
rect test of this hypothesis would be to measure decay rates 
of nuclei placed in a high flux of neutrinos such as produced 
by a nuclear reactor Preliminary work for such an experiment 
has already been carried out at the Pennsylvania State Univer- 
sity Breazeale Nuclear Reactor Although initial results from 
this experiment were inconclusive due to unexpectedly large 
gamma backgrounds, the feasibility of such an experiment was 
established. (It is interesting to note that in 1922 Hevesy [l^ 
searched for modifications of the decay rates of ^^^U and ^"'Pb 
samples that were irradiated with radium; no effect was ob- 
served.) 

The subject of this paper is an alternative approach to search 
for such a new neutrino interaction. A consequence of the 
suggested interaction is that antineutrinos emitted by a sam- 
ple of yS-decaying nuclei should affect the decay rates of the 
remaining nuclei in a sample. Interestingly, the flux of neu- 
trinos that can be achieved in such an experiment (see below), 
is actually greater that the estimated flux in the previously de- 
scribed reactor experiment. In such a case the decay parameter 
A - ln2/ri/2 would become time-dependent, reflecting the de- 
crease in the antineutrino flux accompanying the decrease in 
the number of decaying nuclei. For samples of similar total 
activity, this "self-induced decay" (SID) effect would depend 
on the shape of the sample, since antineutrinos produced in 
different geometries (for example a thin foil a sphere) would 
have different probabilities of influencing the surviving nuclei 
before escaping from the source. The objective of the experi- 
ment described below is to search for such a SID effect, and in 
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the process address the more generic question of whether the 
half-life of source depends on its shape. Although this is obvi- 
ously a fundamental question in ri/2 metrology, it appears not 
to have been specifically addressed previously in the literature 
via a dedicated experiment. Such an experiment could also set 
constraints on other new short-range interactions between nu- 
clei that would affect the decay rates of neighboring nuclei. 

We begin by developing the phenomenology used to charac- 
terize these experiments. According the usual radioactivity de- 
cay law, the activity of a sample containing N(t) unstable nuclei 
of the same type is -dN/dt = -N(t) - AN(t), where /I is a con- 
stant decay rate corresponding to these nuclei. To parametrize 
a possible SID effect, let us assume that the decay parameter A 
for an individual nucleus will depend linearly on the local an- 
tineutrino density Pvir, t). A more general dependence could be 
handled in a similar fashion. We can write the decay rate of a 
nucleus at position ?and time t as A{r, t) - Aq + A\{r, f), where 
Aq represents the intrinsic contribution to the y6-decay rate aris- 
ing from the conventional weak interaction, along with a pos- 
sible constant background arising from new interactions. The 
perturbation A\{f,t) is proportional to the local neutrino den- 
sity, A\(r,t) - ripy(r,t), where the proportionality constant rj 
depends on the strength of the interaction. 

In the presence of the SID effect, the activity of a radioac- 
tive sample no longer obeys the simple decay law; since each 
nucleus experiences a different local neutrino density, the nu- 
clei in the sample do not have the same decay rate. However, 
the extensive experimental work mentioned above shows that if 
the SID effect is real, it must be very small to have gone unde- 
tected, which implies \Ai(r, t)\ <K Aq. This allows one to define 
a sample-averaged decay rate A(t) = Ao + Ai(t), where 



Ai(t) 



Nit) J 



rp{r,t)Ai(r,t). 



(1) 



Here p(r, f) is the number density of undecayed nuclei and the 
integral is over the volume of the sample. To lowest order in rj, 
the neutrino density is proportional to A^(f) so Ai{t) is propor- 
tional to N(t). We can thus write Ai{t) s ^AoNit)/No, where 
N(t = 0) = No, and the dimensionless parameter f depends on 
the geometry of the sample. Notice that ^ - Ai{Q)/Ao repre- 
sents the fractional change in the decay rate at f = due to the 
SID effect. It follows that with Ai{t) the decay law assumes the 
approximate form: 



dN(t) 
dt 



1+^ 



(2) 



Eq. (|2]l can be solved exactly by assuming that A^(f) = 
N{) tx^{-A()t)F{t), in which case Fit) is a solution of the equa- 
tion dF/F^ - exp{-x)dx, where x - AqI. After solving for 
F(x), we find 



N(x)^ 



l+^(l-e-0 



(3) 



We note in passing that Eqs. (|2]i and ^ also describe the non- 
linear growth of a biological population if appropriate signs for 
Ao and ^ are incorporated lll4ll . 



While ^ is a parameter that depends on the geometry of 
the sample under study, we can estimate ^ for a homogeneous 
spherical body of radius R. One can show that for isotropic 
emission, the local density of neutrinos at a distance r from the 
center is 



Pv(r, t) 



Aop 



R R 



(7-^)'"(^) 
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where pit) - 3N{t) / (AnR^) is the average density of undecayed 
nuclei in the sample, and c is the speed of light. (To lowest 
order in rj, we can neglect the spatial variation of p induced by 
the SID effect.) Using this result, we find that the decay rate for 
the spherical sample is given by Eq. (|2]) with ^ — > ^j, where 

For the foil, the corresponding quantity ^/ can be obtained by 
numerical integration. 

We have carried out an experimental search for the SID effect 
by comparing the activities of '^**Au (T1/2 - 2.7 d) as deter- 
mined from a thin Au foil and from a sphere of the same mass 
and similar activity. In this experiment the expression of inter- 
est is the ratio of count-rates for the sphere and foil as a function 
of f = x/Ai). Using Eq. (O, this is given (to first order in ^) by 



A^sphere Nqs 



(1 - -H 2A^e- ') ., 



(6) 



A^foii Nof 

where A^o, and No/ denote the number of '^^Au atoms initially 
present in the sphere and foil, respectively, and A^ = - is 
the difference in ^ for the sphere and foil. A plot of g{x, A^ - 
0.003) is shown for illustration in Fig. [T] where the value of 
A^ has been suggested by the data of Refs. HH H BH]. We 
see from Eq. © that for NoJNof ^ I, g(x = 0, A^) = I + A^, 
while g{x —> 00,^) = 1 - A^. These results can be understood 
by noting that for A^ > and Nos/Nof = 1 the increased decay 
rate arising from the SID effect causes a relative depletion of the 
'''^Au population in the sphere initially compared to the foil, 
such that eventually a crossover in decay rates is achieved at 
t = Ti/2- Thereafter the relatively larger surviving population 
of decaying atoms in the foil produces a larger count rate in the 
foil which persists as f ^ 00. 

As noted above, this experiment can be viewed as a generic 
test of whether nuclear decay rates depend on the shape of the 
sample, irrespective of the underlying mechanism. However, if 
we are specifically interested in whether a shape-dependent SID 
effect could arise from antineutrinos emitted in y6-decay then 
we must account for possible contributions to Aq from the solar 
neutrino flux IF© = 6 x 10'° cm"^s~' . It is evidently desirable to 
start with samples having the largest possible internal antineu- 
trino flux in order to achieve the largest effect. Although the 
initial antineutrino flux in a spherical sample varies as a func- 
tion of the distance r from the center, a useful reference value 
is the flux at the surface r = R which is given by 

\dN/dt\ eApR 



T{R) 



4-kR^ 



5.9xl0^^cm-h-^£R, (7) 



where A - 3. Ox 10 s Ms the '^^Au decay constant, p - 5 



10^'^ cm is the Au atom density, and e = M)('^**Au)/A^()(total). 
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Af = 0.003 



Figure 1 : Plot of g(x, ) from Eq. (6) with = 0.003 and A'o,, /No / = 1 . 

It follows from Eq. O that T^iR) = can be achieved for a 
sphere with R = I mm if e s 1 x 10"'', and these values guided 
the design of our experiment. It should be noted, however, that 
since the effect of neutrinos may depend on their energies as 
well as their flux, choosing parameters such that T^iR) > Tq 
does not necessarily guarantee that the SID neutrinos will give 
the dominant contribution. We note in passing that since the 
specific activity achieved in this experiment exceeds that from 
previous measurements of the '^^Au half-life, the conclusions 
drawn from earlier experiments are not directly relevant in the 
present context. 

Two separate sphere:foil comparisons were performed. 
Squares were cut from pure Au foil (25 fim thick by 1 .4 mm 
square), each of mass 1 mg. Spheres were formed from the 
squares by melting in a borax flux. Each specimen was ir- 
radiated in the RT2 facility [15] of the NIST Center for Neu- 
tron Research for 12 hours at a thermal neutron fluence rate of 
3.1 X 10'^ m^^s"' to produce an activity of 1.1 GBq (30 mCi) 
at the end of irradiation. Under these conditions the contribu- 
tion to the total activity of the double-capture product '^^Au 
(7^1/2 = 3.1 d) is calculated to be 1.5%. For the spherical 1 mg 
sample, 'FiR) = 2.7^0 initially, decreasing to IF© at f = Ti/2- 

The radioactive decay of each sample was followed by re- 
peated counting, using a Ge gamma spectrometry system that 
had been used previously to measure the half-life of '"'^Au El]. 
Sources were held 40 cm from the detector end cap, with an 
absorber of 2.5 cm of Pb between the source and the detector 
Spectra were collected every 2 hr (live time), with 65 to as many 
as 200 spectra in a data set. The initial counting rate was such 
that between 1x10^ and 6x10'' net counts were collected in 
the main gamma peak (41 1.8 keV) in each spectrum. 

Net areas of the three '^^ Au gamma-ray peaks (411 .8, 675.9, 
and 1087.7 keV) were integrated with a fixed-boundary base- 
line subtraction routine fll\. The same peak and baseline chan- 
nels were used for all four data sets. Because of the thick ab- 
sorber, the weaker higher energy lines were enhanced relative 
to the dominant 412 keV line and the low-energy gammas of 
'^^Au were completely undetectable. 

The decay data were fitted to an exponential decay law, cor- 
recting each datum for a linear combination of decay during the 



counting interval, extending dead time (pulse pileup in the ana- 
log circuitry) and non-extending dead time (analog-digital con- 
verter (ADC)) iTs", W]. The fit employed the Solver function 
in Microsoft Excel [20}, a nonlinear reduced gradient method 
for;^'^ minimization to determine the decay constant A, the ini- 
tial activity, and the pileup parameter a. The data points were 
weighted as the inverse square of the Poisson counting preci- 
sion. Individual fits were performed for each of three gamma 
rays in each of four samples (two foils and two spheres). To 
accentuate any effect of the highest specific activity, the twelve 
fits were repeated for a subset of the first thirty spectra, covering 
approximately the first half-life. 

To minimize rate-related instrumental eff'ects, the data were 
first examined in a model-free way, simply comparing peak 
areas of the sphere to the foil in pairs of spectra measured 
at nearly the same counting rate (but slightly diff'erent decay 
times). A small temporal trend in the gamma peak ratios for the 
first experiment was not significantly different from that of the 
associated pulser peak. In the second experiment there was no 
detectable difference in 149 sphere/foil ratios at 412 keV over 
nearly 12 half-lives. 

More fundamentally, half-lives were determined by applying 
the model described in Refs.l 18. T9I1 to each data set: two ex- 
periments, two samples each (sphere vs. foil), and three gamma 
rays per sample. The uncertainty (Is) in the half-life was as- 
signed as the increment that increased the weighted sum of 
squares of normalized residuals by one unit (2I]. The model 
described the data very well: the median per degree of free- 
dom was 1.1 (range 0.9 to 2.6). The results are given in Tablefl] 

To search for effects at the highest specific activity, the decay 
curves were re-fitted using a subset of the data encompassing 
the first half-life. Because of the shorter decay duration, the 
uncertainty of the half-life was greater, as shown in Table |2] 
Half-lives measured in this work are in agreement with those 
previously determined in this laboratory 11611 . 

The measurements in Tables [1] and |2] are dominated by the 
strong 412 keV peak. For the full data sets, the weighted mean 
foil/sphere half-life ratio is 0.9968(14) for experiment 1 and 
0.9999(7) for experiment 2, where the uncertainty is the stan- 
dard deviation of the mean. Assuming a normal distribution, 
the probabilities that the ratio does not differ significantly from 
unity are P = 0.02 and P = 0.90, respectively. The second ex- 
periment, with nearly twice the number of data points, shows 
less difference between the two shapes. 

Our results can be used to estimate by noting from Eq. ^ 
that 



g(x = 0, A^) (A^sphere/A^foil)r=0 



g(x 00, A^) (A^sphere/A^foiOr 



s 1 + 2A^. 



(8) 



For each of experiments 1 and 2 the f = 0(f — > 00) ratio 
was obtained by averaging the 411.8 keV data over five 2- 
hour counts at the beginning (end) of the experiment. We find 
for the ratio in Eq. ^ the value 0.9964(7) from experiment 
1, and 0.9992(16) from experiment 2; the weighted mean of 
both experiments gives 0.9969(10), where the Icr uncertain- 
ties are obtained from propagated Poisson statistics. Insert- 
ing these results into Eq. ^ yields A^i - -1.8(4) x 10"^, 
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Table 1: Half-lives (in hours) determined for each gamma-ray in each sample, using all n data points. The symbol ± here and throughout this work denotes the 
standard (Is) uncertainty. 



Energy (keV) 


Foil Fl (« = 99) 


Sphere SI (n = 65) 


Foil F2 (n = 146) 


Sphere S2 (« = 197) 


412 


64.517 + 0.019 


64.709 + 0.046 


64.573 + 0.016 


64.589 + 0.009 


676 


64.504 + 0.061 


64.895 + 0.150 


64.554 ± 0.040 


64.562 + 0.024 


1088 


64.705 + 0.095 


64.790 ± 0.220 


64.716 + 0.057 


64.584 + 0.033 


7-Foil /T^Sphere 
-' 1/2 ' -' 1/2 


0.9968 + 0.0014 


0.9999 + 0.0007 



Table 2: Half-lives (in hours) determined from the initial n = 30 spectra (first 63 hours). 



Energy (keV) 


Foil Fl 


Sphere SI 


Foil F2 


Sphere S2 


412 


64.92 + 0.24 


64.38 + 0.24 


63.69 + 0.39 


64.82 + 0.44 


676 


63.76 + 0.73 


64.12 + 0.73 


63.77+ 1.05 


66.82+ 1.24 


1088 


64.32 + 1.08 


65.08 + 1.20 


65.12+1.49 


66.84+ 1.73 


T-Foil /T^Sphere 
-' 1/2 ' -' 1/2 




1.006 + 0.008 


0.979 + 0.008 





A^2 = -4(8) X 10""*, and A^,„ = -1.55(50) x 10"^ obtained 
from experiments 1, 2, and their mean, respectively. For pur- 
poses of comparison, the peak-to-trough variation in the BNL 
Si/Cl ratio of count rates, Ar fst] is ~ 3 x 10"^, corresponding 
to a change A/ in solar flux of 6.683 x 10"^ from perihelion to 
aphelion. The ratio Ar/A/ w 4x10"^ provides a measure of the 
improved sensitivity of the current experiment when compared 
to the above values of A^. 

Although the preceding comparisons of the f = and f — » oo 
count-rates exhibit small deviations from the expected null re- 
sults, the overall half-lives measured in this work are in agree- 
ment with those previously measured in this laboratory [16]. 
From the full data sets the mean foil/sphere half-life ratio is 
0.999 + 0.002, and for the initial decay period the ratio is 
1.00 + 0.02. These results are consistent with the deviations 
previously noted in A^i and A^,„: returning to Fig. [T]we note 
that there is a degree of cancellation that arises when T[/2 for 
the spherical samples is obtained by combining data for t <T\i2 
with data for t > Ti/2- Separating the data from these two time 
intervals is thus always desirable, but comes at a cost in statis- 
tics, as we have noted previously. 

In summary, the present experiment is the first direct preci- 
sion test of whether the decay rate of a radioactive source de- 
pends on its shape. Our results in Table [T]indicate a ~2.3cr de- 
viation of the foil/sphere ratio in experiment 1 from unity. From 
Table|2] based on the initial 30 spectra, the foil/sphere ratio for 
experiment 2 deviates from unity by ~2.6cr. These results thus 
leave open the possibility that the half-life of a radioactive nu- 
clide could in fact depend on its shape (due to the internal flux 
of neutrinos, photons, or electrons), and hence suggests that 
additional experiments are necessary. We are thus planning to 
repeat this experiment in an attempt to reduce our uncertainties 
by starting with foil and sphere samples with greater activity. If 
a significant eff'ect attributable to neutrinos were seen in such 
an experiment, this would support the inference drawn from 
Ref. SB 0,11101 of a possible solar influence on nuclear decay 
rates. 
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